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Abstract 

The aim of this article is to demonstrate how the vector field method of Klain- 
erman can be adapted to the study of transport equations. After an illustration 
of the method for the free transport operator, we apply the vector field method 
to the Vlasov-Poisson system in dimension 3 or greater. The main results are 
optimal decay estimates and the propagation of global bounds for commuted 
fields associated with the conservation laws of the free transport operators, un¬ 
der some smallness assumption. Similar decay estimates had been obtained 
previously by Hwang, Rendall and Velazquez using the method of characteris¬ 
tics, but the results presented here are the first to contain the global bounds 
for commuted fields and the optimal spatial decay estimates. In dimension 4 
or greater, it suffices to use the standard vector fields commuting with the free 
transport operator while in dimension 3, the rate of decay is such that these 
vector fields would generate a logarithmic loss. Instead, we construct modified 
vector fields where the modification depends on the solution itself. 

The methods of this paper, being based on commutation vector fields and 
conservation laws, are applicable in principle to a wide range of systems, in¬ 
cluding the Einstein-Vlasov and the Vlasov-Nordstrom system. 
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1 Introduction 

A standard approach to the study of asymptotic stability of stationary solutions of 
non-linear evolution equations consists in an appropriate linearization of the sys¬ 
tem togetheQ with 

1. a robust method for proving decay of solutions to the linearized equations, 

2. an appropriate set of estimates for the non-linear terms of the original system, 
using the linear decay estimates obtained previously. 

For systems of non-linear wave equations such as the Einstein vacuum equations 
Ric{g ) = 0, several methods for proving decay of solutions to the linear wave equa¬ 
tion □)// = 0 where □ = -d 1 2 3 t + ^ d 2 , is the wave operator of the flat Minkowski spac^f] 

i=i 

are a priori available. One of the classical methods to derive decay estimates is to use 
an explicit representation of the solutions, such as the Fourier representation, to¬ 
gether with specific estimates for singular or oscillatory integrals. While this method 
provides very precise estimates on the solutions, it does not seem sufficiently robust 
to be applicable to quasilinear system of wave equations such as the Einstein equa¬ 
tions, and the method of choicqjfor proving decay in view of such applications is the 
commutation vector field method of Klainerman (7] and its extensions using mul¬ 
tiplier vector fields, see for instance 11311141 [3l . The method of Klainerman is based 
on 

1. A coercive conservation law: the standard energy estimate in the case of the 
wave equation. 

2. Commutation vector fields: these are typically associated with the symmetries 
of the equations. In the case of the wave equation, these are the Killing and 
conformal Killing fields of the Minkowski space. 

1 A third ingredient not needed in the present case is that of modulation theory, see for instance [8] for 
an application of modulation theory in the context of the Vlasov-Poisson system. 

2 In case of perbutations around a non-flat solution with metric g, the operator □ would naturally be 
replaced by HHg, the wave operator of the metric g. 

3 More recently, a mix of microlocal and vector field methods have also been successfully developped, 
in particular to handle complex geometries involving trapped trajectories, see for instance [12] for an 
application of these tools. 
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3. Weighted vector field idendities and weighted Sobolev inequalities: the usual 
vector fields d t ,d x i are rewritten in terms of the commutation vector fields. 
The coefficients involved in these decompositions contain weights in t and 
\x\ and the presence of these weights leads to weighted Sobolev inequalities, 
that is to say decay estimates. 

The typical method used in the study of the Vlasov-Poisson and other systems 
of transport equations such as the Vlasov-Nordstrom system is the method of char¬ 
acteristics. This is an explicit representation of the solutions and thus, in our opin¬ 
ion, should be compared with the Fourier representation for solutions of the wave 
equation. What would then be the analogue of the vector field method for trans¬ 
port equations? The aim of this article is twofold. First, we will provide a vector 
field method for the free transport operator. In fact, in a joint work with J. Joudioux 
and D. Fajman, we have developped a vector field approach to decay of averages 
not only for the free (non-relativistic) transport operator but also for the massive 
and massless relativistic transport operators, see (4j. In this paper, we will give two 
different proofs of Klainerman-Sobolev inequalities. The easier proof will give us 
a decay estimate for velocity averages of sufficiently regular distribution functions, 
i.e. quantities such as / veH „ f{t,x, v)clv. However, this proof fails in the case of ve¬ 
locity averages of absolute values of distribution functions, i.e. quantities such as 
IveR n I f\(.t,x, v)dv, because higher derivatives of |/| will typically not lie in L 1 even 
if / is in some high regularity Sobolev space. On the other hand, the decay estimate 
obtained via the method of characteristics can be applied equally well to / and |/|. 
We shall therefore give a second proof of Klainerman-Sobolev inequalities for veloc¬ 
ity averages which will be applicable to absolute values of regular distribution func¬ 
tions. The first approach, which is closer to the standard proof of the Klainerman- 
Sobolev inequality for wave equations, consists essentially of two steps, a weighted 
Sobolev type inequality for functions in L\. and an application of this inequality to 
velocity averages, exploiting the commutation vector fields. The improvement in 
the second approach comes from mixing the two steps together. 

In the second part of this paper, we will apply our method to the Vlasov-Poisson 
system in dimension n > 3 


d t f+v.\7 x f + pV x (p.V v f = 0, 

(1) 

A (p = p{f), 

(2) 

Of 

II 

o 

II 

o” 

(3) 

n 

where p = +1, A = - ^ d 2 xi , f = f{t,x, v) with t £ IK, x, v e IR", 

/o is a sufficiently 


i= 1 


regular function of x, v and p{f) is given by 

p(f){t,x):~ f f(t,x,v)d n v. 

Jve R" 

Our main result can be summarized as follows (a more precise version is given 
in Section l4~2l . 

Theorem 1.1. Let n>3 and N > ^ + 2 ifn > 4 and N > 14 ifn — 3. Let 0 <5 < 

Then, there exists cq > 0 such that for all 0 < e < eo, ifEfjg[fo] ^ e, where E/v.gl/ol is 
a nor /fQ containing up to N derivatives offo, then the classical solution f{t,x, v) of 

4 See Section [4711 for a precise definition of the norms. The 8 encodes some additional integrability 
properties of the solutions. 
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GMS exists globally in and satisfies the estimates, V t e R andMx e l 

1. Global bounds 

E NiS [f](t)<2e. 

2. Space and time pointwise decay of averages of p(f) 

for any multi-index a with |a| < TV- n, 

Cn,ji,s £ 


( 4 ) 


\p(Z a f)(t,x )\< 


(1 +|f| +|JC|)" 


where Z a is a differential operator of order a obtained as a combination of\a\ 
commuting vector fields and Cm,„,s > 0 is a constant depending only on N,n,S. 


3. Improved decay estimates for derivatives of / 

for any multi-index a with |a| < TV- n, 


Cjv,n,<5 e 


(l + |r| + |x|)” +l “ r 
4. Boundedness of the L l+S norms of V 2 </> and V 2 Z“</> 

for any multi-index a with \a\ < TV, [|V 2 Z“0(t)|| i i +< 5 {H n) < CN,n,S e - 


5. Space and time decay of the gradient of the potential and its derivatives 
for any multi-index a with \a\ < AT- (3n/2 + 1), 

\VZ a f{t,x )\< ■ 

as well as the improved decay estimates 

Cn,ii,8E 


CN,n,fi£ 

tfr-W2(l + \t\ + \ X \) nl2 ' 


|d“V^)(t,x)| < 


t^-2)l2 [l + \ t \ + \ x][) nl 2 + \a\- 


Remark 1.1. Stronger bounds can be propagated by the equations provided the data 
enjoy additional integrability conditions. More precisley, the improved decay esti¬ 
mates for derivatives of p{f) can be improved to 


\p(d\d a x f){t,x)\ < 


_Q.yC_ 

(l + |r| + |x|)" +l “ l+T 


and the improved decay estimates for derivatives of the gradient off can be improved 
to 


\d r t d a x \/f(t,x)\< 


_Cjye_ 

ft"-2)/2( 1 + |r| + |x|)”/2+l“l+T’ 


the point being that additional t derivatives now bring additional decay in t and\x\. 
These stronger estimates hold provided the initial data have stronger decay in x, v 
than what is needed in the proof of Theorem li.il Similarly, one can propagates If 
norms with p > 2 for V<p andV Z a (j) provided additional v decay of the initial data is 
assumed. 

5 Under some mild conditions on the initial data, global existence is already guaranteed from the works 
|T7j[II] , so the main points of the theorem, apart from providing an illustration of our new method, are 
the propagation of the global bounds and the optimal space and time decay estimates for the solutions. 
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Remark 1.2. Similar time decay estimates have been obtained in ©/or derivatives 
ofpif ) and <p using the method of characteristics under different assumptions on the 
initial data. On the other hand, the optimal decay rates in space and the propagation 
of the global bounds (4) were, as far as we know, not known prior to our work. 

Remark 1.3. As' is clear from the proof below and is typical of strategies based on com¬ 
mutation formulae and conservation laws, the method is very robust. In particular, 
we are not using the method of characteristics, nor the conservation of the total energy 
for the system (T)-(3) . An illustration of this robustness will be given in © where we 
will apply a similar approach to the study of the Vlasov-Nordstrom system. 

Previous work on the Vlasov-Poisson system and discussion 

There exists a large litterature on the Vlasov-Poisson system. We refer to the intro¬ 
duction in [L5]| for a good introduction to the subject and only quote here the most 
important results from the point of view of this article. In the pioneered work (TJ, 
small data global existence in dimension 3 for the Vlasov-Poisson system was es¬ 
tablished together with optimal time decay rates for p{f) and V/ but no decay was 
obtained for their derivatives. The optimal time (but not spatial) decay rates for 
derivatives of p(f) and V(p has been only much later obtained in (6), covering at the 
same time all dimensions n > 3. Both these works use decay estimates obtained 
via the method of characteristics. In fact, in [1] and even more in (6) , precise es¬ 
timates on the deviation of the characteristics from the characteristics of the free 
transport operator are needed in order to obtain the desired decay estimates. Par- 
allely to these works giving information on the asymptotics of small data solutions, 
let us mention that under fairly weak assumption on the initial data (in particular, 
no smallness assumption is needed), it is known that global existence holds in di¬ 
mension 3 for the solutions of G)-(3) , see nzjnn- The strongest results concerning 
the stability of non-trivial stationnary solutions of (TJ - (3) with p = -1 have been 
obtained in (8). They are not based on decay estimates but on a variational charac¬ 
terisation of the stationary solutions. On the other hand, this type of method does 
not provide asymptotic stability of the solutions but orbital stability. It is likely that 
any result addressing the question of asymptotic stability will need to go back to an 
appropriate linearization of the equations combined with robust decay estimate^]. 
We believe that, once again, the vector field method would be totally appropriate for 
the derivation of such decay estimates. Finally, let us mention the celebrated work 
IM on Landau damping concerning the stability of stationnary solutions to {TJ - (3} 
with periodic initial data. In view of the present work, it will be interesting to try to 
revisit this question using vector field methods. 

Outline of the paper 

In Section[2] we introduce the vector fields commuting with the free transport op¬ 
erator and the notations that we will use throughout the paper. In Section [3J we 
present and prove decay estimates for velocity averages of solutions to the transport 
equation. In the following section, we present our results on the Vlasov-Poisson sys¬ 
tem. The remaining last two sections are devoted to the proof of these results, first 
in dimension n > 4 and then in dimension 3 using modified vector fields. 

6 See for instance [5] for some stability results using the linearization approach in the case of the 
spherically-symmetric King model. 
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2 Preliminaries 

Throughout this article, / will denote a sufficiently regular function of ( t, x, v ) with 
t e [R and (x, v) eR" x M”. By sufficiently regular, we essentially mean that / is such 
that all the terms appearing in the equations make sense as distributions and that 
all the norms appearing in the estimates are finite. For simplicity, the reader might 
just assume that / is smooth with compact support in x, v (but any sufficient fall-off 
will be enough). 

We will denote by T the free transport operator i.e. 

T{f):=d t f+f i v * i d x ,f, 

i =1 

where d t f = and for all 1 < i<n, d x if = Similarly, for any sufficiently regular 
scalar function (p , l',], will denote the perturbed transport operator 

r 0 (/):=r(/)+pv^.v y /, (5) 

where ju = ±1, corresponding to an attractive or repulsive force. Since we are dealing 
only with small data solutions, the sign of /i will play no role in the rest of this article. 

The notation A< B will be used to specify that there exists a universal constant 
C > 0 such that A < CB, where typically C will depends only on the number of di¬ 
mensions n and a few other fixed constants, such as the maximum number of com¬ 
mutations. 

2.1 Macroscopic and microscopic vector fields 

Consider first the following set of vector fields 

• Translations in space and time d t ,d x i, 

• Uniform motion in one spatial directioif] td x i , 

• Rotations x l d x j-x ] d x i, 

n 

• Scaling in space ^ x l d x i , 

i =1 

7 Recall that these vector fields are the generators of the Gallilean transformations of the form x e U n — 

x+ t Vi, where v f* = 8 

l ’ 11 
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n 

• Scaling in space and time td t + y x 1 d x i. 

i =1 

The above set of vector fields is associated with the Gallilean invariance of macro¬ 
scopic fields and equations. We will denote by f the set of all such vector fields 

n n 

r = \ d t ,d x i, td x i,x l d x j -X ] d x i, £ x'd xi ,tdt + J^ x l d x i, 1 <i,j<n\. 

I (=1 i=i J 

One easily check that while the translations commute with T, uniform motions, 
rotations or the scaling in space do not. The correct replacement for these vector 
fields is most easily explained using the language of differential geometry; the inter¬ 
ested reader may consult [18] 0] for detailed constructions (in the case of the rela¬ 
tivistic transport operator). We shall here only present the resulting objects which 
are the vector fields 


• Uniform motions in one direction in microscopic form td x i +d v i, 

• Rotations in microscopic form x l d x j - x J d x , + v l d v j - v 1 d v i, 

n 

• Scaling in space in microscopic form ^ x l d x i + v'd v i. 

i=i 

One can then easily check 

Lemma 2.1 (Commutation with the transport operator). 

• If Z is any of the translations, microscopic uniform motions or microscopic ro¬ 
tations, then [T,Z] - 0. 

• IfZ is the microscopic scaling in space, then [T,Z] = 0. 

• IfZ is the scaling in space and time, then [T,Z] = T. 

Remark 2.1. From the two scaling commuting vector fields, it follows automatically 
that td t - X” =1 v l d v i also commutes with T in the sense that 

[T,td t -Y j v i d v i] = T. 

i=i 

This vector field will be used to obtain improved decay for t derivatives of velocity 
averages. 

To ease the notation, we will denote by OT := x‘ d x j - x J d x , the rotation vec¬ 
tor fields in x and by Q v .. v l d v j - v 1 dy the rotation vector fields in v. The full 
microscopic rotation vector fields are thus of the form Qf. + Of... Similarly, we will 

n 

denote by S x + S v := ^ x l d x i + v l d v i the scaling in space in microscopic form, with 

i=i 

n n 

S x = y x'd x i and S v = ^ v'd v i. 

i=i i=i 

Let now y be the set of all the above microscopic vector fields including the trans¬ 
lations and the space and time scaling i.e. 


j=\d t ,d x i,td x 


,n*. + n y .,s x 


- S v , td t + ^ x l d x i, 1 < i,j < n 

i=i 
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2.2 Multi-index notations 


Let Z l , i = 1,..,2n + 3 + n[n - l)/2 be an ordering of L. For any multi-index a, we 
will denote by Z a the differential operator of order |a| given by the composition 
Z ai Z“ 2 ... 

In view of the above discussion, to any vector field of F, we can associate a unique 
vector field of y. More precisely, 

d t - d t . 


d x‘ 

- d x‘> 

td x‘ 

- td x i +d v i, 

nf. 

— QL+nf,, 

i t 

l J l J 

s x 

- s* + s v , 

n 

n 

td t + Y,x l d x i 

—|• td t + ’^x l d x , 


i= 1 i'=l 

Thus, to any ordering of F we can associate an ordering of y. We will by a small abuse 
of notation, denote again by Z l the elements of such an ordering since it will be clear 
that, if Z l is applied to a macroscopic quantity, such as a velocity average, then Z 1 e 
T and if Z ! is applied to a microscopic quantity, i.e. any function depending on (x, v) 
(and possibly t), then Z 1 e y. Similarly, for any multi-index a, we will also denote by 
Z a the differential operator of order \ a\ given by the composition Z“’ Z“ 2 .. obtained 
from the vector fields of y. 

For some of the estimates below, it will be sufficient to only consider a subset 

of all the vector fields of T and y. Let us thus denote by F s the set of all the macro- 

n 

scopic vector fields apart from d t and td, + ^ x l d x i, which are the only vector fields 

i= 1 

containing time derivatives and by y s the corresponding set of microscopic vector 
fields, i.e. 


r 5 = | d x i , td x i , x l d xj - x J d x i,Y,x l d x i, 1 < i,j < n 

T = {d x i,td x i+d v i,Q. x ij +n v ij ,S x + S v ,\<i,j<ny 

The notation Z“ e r'“ (respectively Z“ e y 1 * 1 ) will be used to denote a generic 
differential operator of order |a| obtained as a composition of \a\ vector fields in F s 
(respectively in y s ). The standard notation d x will also be used to denote a differen¬ 
tial operator of order \a\ obtained as a composition of \a\ translations among the d x , 
vector fields. 

The following lemmae can easily be checked. 

Lemma 2.2 (Commutation within F). For any Z a e y'“', Z a ' e y 1 "' , where a, a' are 
multindices, we have 

[Z“,Z a ']= £ c a ’ a 'zP, 

|/3|<|a|+|a'|-l P 

for some constant coefficients c“’“ . Moreover, if Z a ,Z a ' e y s , then all the ZJ’’ of the 
right-hand side belongs to y s . 
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Lemma 2.3 (Commutation of Z“ and weights in v). Let q > 0. For any sufficienly 
regular function f of{ t, x, v) and for any Z a e where a is a multi-index, we have 

|Z“ [(1 + v 2 ) ql2 f] | < (1 + v 2 ) ql2 £ |z^(/)|. 

\P\S\a\ 

Moreover, ifZ a e then all the Z& belong to y|f 1 in the above inequality. 


2.3 Velocity averages and commutators 

For any integrable function / of v e M'\ we will denote by p(/) the quantity 

p(f):=[ fd n v 
JveR n 

We have the following lemma. 

Lemma 2.4. For any sufficiently regular function f of(x, v ) we have 

• for all 1 < i < n, 

d xi p(f)=p(d x if), 

• for all t eU and all l < i < n, 

td x ipif) = p([td x i +d v i)f), 


• for all 1 < i,j < n, 


n 


= p ((nfj+(/)). 


for allteU and tel", 


i =1 


tt 


td t + £ x'd xi p(/) = p + £ x ! d xi (/) 


VV i=t 


• andfinally 

S x p{f) = p[[S x + S v ) lf)) + np(f), 

where S x and S x + S v are the spatial scaling vector fields in macroscopic and 
microscopic forms. 


Proof. The proof is straigtforward and consist in identifying total derivatives in v. 
For instance, we have 


>(s y )(/)= f 

Jv 


, L v ‘ d v‘f 

veR" \ i = i 


d n v- f -nfd n v, 

JveR n 


where we have integrated by parts in each of the v'. Similarly, in the case of rota¬ 
tions, it suffices to note that for any 1 < i, j < n, QC is an angular derivative in v and 
thefore, f veK n Fl v .. [f)d n v = 0. □ 

In the remainder of this paper, we shall write the preceding lemma as 


Zp(/) = p(Z(/)) + c z p{f), 

where we are using, by a small abuse of notation, the letter Z to denote a generic 
macroscopic vector field and its corresponding microscopic version and where cz = 
0 unless Z is the spatial scaling vector field, in which case cz - n. 
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2.4 Vector field identities 

The following well-known identity will be used later 
Lemma 2.5. For any 1 < j < n, we have 


\x\ 2 d xj = Y J X i nF+x J S x , 

i= 1 


( 6 ) 


and thus, at any x + 0, 


n r 1 vl 

I x\d xJ = Y—nf i + —s x 
Xl tl 1*1 1 1*1 


where the coefficients are all homogeneous of degree 0 and therefore uniformly 
bounded. 


The following higher order version will be used often in the derivation of the 
Klainerman-Sobolev inequalities of the next section. 

Lemma 2.6. For any multi-index a, 

(f+ |x|)“d“ = £ CpZ?, 

\p\<\a\,zP£rf' 

where the coefficients Cp are all uniformly bounded. 

Proof. The lemma is a consequence the previous decomposition, the fact td x i is 
part of our algebra of commuted vector fields and that \d x i ,t+\x\ is homogeneous 
of degree 0. □ 


2.5 The commuted equations 

We now turn to the study of the transport operator T,p defined by (5). Many of the es¬ 
timates below are only valid provided <p has sufficient regularity. In the applications 
to the Vlasov-Poisson system of this article, we will eventually control the regularity 
of <p via a bootstrap argument. For all the estimates below, we therefore assume that 
cf> is a sufficiently regulaiQ function of ( t, x) defined on [0, T] xij, for some T > 0, 
which decays sufficiently fast as \x\ —► +oo. 

The following lemma can then easily be checked. 

Lemma 2.7. Let f be a sufficiently regular function of (t,x,v) and let a be a multi- 
index. Then, there exists constant coefficients such that, 

[T<p,Z a ]f= £ £ Cg V x Z r (p.SJ V Z^f. 

|/3|<|a|-l, |y|+|/5|<|a| 

Moreover ifZ a e T^', then Z r e T^ and Z@ e y!f in the above decomposition. 
Similarly one has 

8 For instance, one can assume that (p is a smooth function on [0, T] x Rj with compact support in x. 
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Lemma 2.8. Letf be a sufficiently regular function of{ t, x, v ) and lety/ be the solution 
to the Poisson equation Ai// = p(/). Then, for any multi-index a, Z a y/ is solution to 
an equation of the form 

AZ>= £ C“Z^p(/), (7) 

ipi<|a| 

where C'* are constants. 

Proof. The lemma is an easy consequence of the fact that all macroscopic vector 
fields apart from the two scalings commute with A while for the spatial scaling S x 
and the space-time scaling td t + S x we have [A, S x ] = 2A and [A, td, + S*] = 2A. 

□ 


2.6 Conservation laws 

We shall use the following (approximate) conservation laws. 

Lemma 2.9. For any sufficiently regular function f of {t,x, v), we have, for all t e 
[0,71, 

ll/(0lll, 1 (R'’xR'|) - ll/(0)lll 1 (RjxR' i ;) \\T<l)(.f){s)\\ L l^n xR n- > ds. 

Similarly, we have for all p > l, for all t e [0, T], 

- H/ (0 ) IIip(r;xC) + Pl ll/ p_1 

and for all q > 1, 

||(1 + V 2 ) C,l2 f P {t)\\ L l( R n xR n ] < ||(1 + V 2 ) ql2 f P {0)\\ L l {R n xR n ] 

+ J 0 11(1+ l. 2 ) /2 /P-! (/) (5) 11il (R ? X R«, 


+ [ ||/ p (l + v 2 ) [q 1 ), 2 d x (p(.s)\\ LHR n xK) ds. 

•J 0 


( 8 ) 


Note in particular that the conclusions of the lemma hold true when T,/, - T (i.e. when 

d x <p = o;. 

Proof. These are classical estimates so we only sketch their proofs. 

One has (in the sense of distribution), 

[(1 + v 2 ) ql 2 \f\ p ] < (1 + 7) p/2 |/l p_1 |T 0 (/)| + (1 + v 2 ) [q ~ l)l2 \d x <p\\f\ p . 

Using a standard procedural, one can regularize the previous inequality. We shall 
therefore neglect regularity issues here. Integrating the previous line in (f, x, v ) leads 
to 

m r 0 [(l+ v 2 ) ql 2 \f\ p ] dvdxds < 

' 

[(1 + v 2 ) ql2 \f\ p ~ l \T ( p{f)\ + (1 + v 2 ) iq - 1V2 \d x (p\\f\ p ] dvdxds. (9) 


/ 7/1 


9 For instance, assume first that f has compact support in (x, v) with a uniform bound on the support 

of / in (x, v ) for t e [0, T]. For all e > 0, consider the function f e = yje 2 + / 2 ^(x, v), where j(x, v ) is a 
smooth cut-off function which is 1 on the support of / and vanishes for large x and v. Apply then the 
previous estimates to f £ and take the limit e —*► 0. A standard density argument deals with the case of 
non-compact support. 
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On the other-hand, remembering that 7)/, = d t +Y _? i v’d l x +[iV x (p.\7 v and integrating 
by parts in x and v, we obtain 



||(1 + y 2 )'7 /2 /P(t)|| i i (R „ xR „ ) 


which combined with (9} leads to the desired estimate (8). 


□ 


3 Decay of velocity averages for the free transport oper¬ 
ator via the vector field method 

Since the main purpose of this article is to illustrate how the vector field method can 
lead to robust decay estimates for velocity averages, let us for the sake of comparison 
recall the Bardos-Degond decay estimate and its proof. 

Proposition 3.1 (Q]). Let f be a sufficiently regular solution ofT(f) = 0. Then, we 
have the estimate, for all t > 0 and all rel", 

\p{f)\{t,x) < f sup |/(0, x,v)\d n x. (10) 

t JxeR n deR" 

Proof. The proof of this classical estimate is based on the method of charateristics. 
More precisly, if / is a regular solution to T{f) = 0, then it follows that 


f{t,x, v) = /(0, x- vt, v ) 


for all teR, x,veR n . 
We then have 



Applying now the change of coordinates y = vt for t > 0 leads to 



□ 


In the above proof, the two key ingredients are 

• the explicit representation obtained via the method of characteristics, 

• the change of variables y-vt. 
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Note that in the presence of a perturbation of the free transport operator, to exploit a 
similar change of variables would require estimates on the Jacobian associated with 
the differential of the characteristic flow, see mm. 

Let us now show how the vector field method can be used as an alternative to 
obtain similar decay estimates. As explained in the introduction, we will give two 
different proofs. 

The first proof will give us decay estimates for quantities of the form f v£R n fd v. 
The starting point of this approach is the following Klainerman-Sobolev inequality 
using L 1 (R"J norms of commuted fields. 


Lemma 3.1 (L ] Klainerman-Sobolev inequality). For any sufficiently regular func¬ 
tion y/ defined on IR", we have 


WW < 


1 

(l + t+\x\) n 


Y ||Z“(Vr)[| L i (H!) . 

\a\<n,Z a eT s 


Proof. This is relatively standard material and we adapt here the presentation given 
in (T9) Chap.2 to our setting. 

Fix (i,x)eR ( x IR” and let y/ be the function 

y/:U n — U (11) 

y —■ yKy):=yr(x+(t+\x\)y). (12) 

Applying a standarc0 Sobolev inequality, we have 

\y/(t,x)\ - |^(0)] < Y I 1 5 “^Hi 1 (b k ( 0 ,i/ 2 )). (13) 

\a\<n 

where B n {0, 1/2) denote the ball in R'! of radius 1/2. On the other hand, we have 
dyiy/ly) = it + \x\)d x iy/(x+ (f + |x|)y) and thus, for y e B n {0,l/2), 

\d“y/[y)\ < Y (.t+\x\) l3 \d^y/[x+{t+\x\)y)\, 

W\s\a\ 

< Y {t+\x+{t+\x\)y\)P\dPy/[x+{t+\x\)y)l, 
l/3|<|a| 

< Y \zPy/(x + {t + \x\)y)\, 

\P\<\a\,Zf>ET s 

where we have used Lemma f2T6l in the last step and the fact that t + \x\ and t + 
\x+ (t + |x|)y| are comparable for y e B n {0, 1/2). Inserting the last line in the Sobolev 
inequality fl3l and applying the change ofvariables z = (f + x|)_y conclude the proof 
of the lemma. 

□ 


Using Lemma [2741 we now note that for any vector field Z and any sufficiently 
regular function / of x, v, we have 


II^(P(/))IIu(r;) = Wp(Zlf)) + czp{f)\\ L i m n x) 

and thus, 

II-^(p(/))IIl i (r;) ^ II^(/)IIli(r;xrs + llp(/)llz. 1 ( [R " x|R ”)- 

Combined with the previous Klainerman-Sobolev inequality, we obtain 

10 Recall that, while the general W k, P( R”) •—> L°°(R n ) embedding requires k > the special case p = 1 

only needs k> n. 
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Proposition 3.2 (Global Klainerman-Sobolev inequality for velocity averages). For 
any sufficiently regular function f defined on 05” x 05”, we have, for all t > 0 and all 
x e 05”, 

lp(/)l(x)r$ / £ \\Z a f\\ LHRM (14) 

U + f+|X|) \a\<n,Z«ey^ 

Note that the above inequality cannot be apply to |/| even is / is say a smooth, 
compactly supported function. Indeed, if |/| is say in W n,p , then |/| is in W 1,p but, 
unless / has some extra special properties, |/| € W n,p when n > 2. On the other 
hand, dTOt clearly holds both for / and for |/|. 

Two disctinct steps lead to the proof of (M) , the L 1 Klainerman-Sobolev inequal¬ 
ity of Lemma l3.1l and the special commutation properties of the velocity averaging 
operator as described in Lemma 12.41 To improve upon JT~T1 , the strategy is to try 
to use at the same time arguments similar to those of Lemma 13.11 and Lemma l2~4l 
instead of applying them one after the other. This will lead to us to the following 
improvement. 

Proposition 3.3 (Global Klainerman-Sobolev inequality for velocity averages of ab¬ 
solute values). For any sufficiently regular function f defined on 15” x IJ5”, we have, 
for all t> 0 and all x e 15”, 

lp(l/l)IM< 1 „ £ \\Z a f\\ L HuM (15) 

(l + f+'x!) |a|snZaer N 

Proof. Let us assume that / is smooth and compactly supported for simplicity. 

Define yj similarly to fTTI as 


y/:B n (0,1/2) 

y 




:= f \f\(x+(t+\x\)y,v)d n v. 
Jve R n 


Next, recall that for any y/ e W 1,1 , \y/\ e W 1,1 and d\y/\ = j^dy/ (in the sense of distri¬ 
bution, see for instance (9) , Chap 6.17) so that in particular \d\y /\| < \dyr\. Let us write 
y = (y| ,..,}'„) and let 8 - Using a 1 dimensional Sobolev inequality, we have 

Hp(0)|< f {\d yi y/{y 1 ,0,..,0)\ + \y}(y 1 ,0,..,0)\)dy 1 . (16) 

dlvi |<<5 1/2 


Now, 


<9 yi ip(yi,0,..,0) = f (t+ |x|)d Xl |/| (x+ (t+ |x|)(yi,0, ..,0)) d n v, 

Jve H" 

= f td Xl \f\ [t,x+ [t + |x|)(yi,0,..,0)] d n v 

Jve H n 

+ f \x\d Xl |/| (x+ (t+ |x|)(yi,0,..,0)] d n v. 

Jvett" 

As before, we can introduce total derivatives in v. For instance, 

f td Xl |/| (x + (f + |x|) (yi, 0,.., 0)) d n v = 

Jve R" 

f [td Xl +d Vl )\f\[x+{t+\x\)(y 1 ,0,..,0))d n v, 

JveU n 

< f \[td Xl +d v f) f[t,x+{t+ |x|)(yj,0,..,0))| d n v, 

Jve R« 


(17) 
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and, using that (yi,0,..,0) e B n { 0,1/2) if |yi| < S 112 , a similar argument can be used 
to handle the second term on the right-hand side of IT71 . This gives us 

|5 yi t//(yi,0,..,0)| < Y f \Zf[x+{t+ |x|)(yi,0,..,0)) | d n v. 

z s E rs Jvm n 

Combined with fl6t . we have obtained that 

IV(0)l $ V f f \Zf[ x + (f + |x|)(yi,0,..,0)) | d n vdy\ 

Z s £Y s J\yi\sS 112 JvzR" 

+ [ f |/(x+ (t + |x|)(yi,0, ..,0)) I d n vdy\. 

J\y\\<S 1 ' 2 J weR" 

We now repeat the same argument varying the variable y 2 - 


|i/d0)|<( f 

J\yi\<8 112 J|y 2 |<5 1/2 

Y [ \Z a f[x+(t+\x\){yi,y 2 , 0 ,.., 0 ))\d n vdy 1 dy 2 . 

|a|<2, JveU” 

Z a zy[ a[ 


Iterating again this argument until all variables y,- appears in the integral of the right- 
hand side, we obtain 


l^(0)|< 


l 


Y [ \Z a f[x+(t+\x\)y)\d n vdy } 
ye-B«(0,1/2) \a\<n, Jve.M n 


Z a £ T s 


a | 


and the conclusion of the proof follows as in the proof of Lemma [3Tl bv the change 
ofvariable z - (t + |x|)y. □ 


We now recall that if / is a solution to T (/) = 0, then, in view of the commutation 
properties of Lemma 12.11 and standard properties of differentiation of the absolute 
value, so are |/| and the commuted fields |Z“/|. Thus, from Lemma 12.91 all the 
norms on the right-hand side of fi~5l are preserved by the flow and we have obtained 
the decay estimate 

Proposition 3.4 (Decay estimates for velocity averages). For any sufficiently regular 
solution f to T(f) - 0, we have, for all t > 0 and all tel", 

\pi\f\)\(t,x)<, E \\Z a nt = m LHKM (18) 

U + t+|X|) |8| <n,Z«^f 

Remark 3.1. If we restrict the set of vector fields only to the uniform motions td x i +d v i, 
then we still obtain the time decay estimate, 

\p{\f\m,x)<^ Y ii^7(f = o)ii L i Mi;) . (i9) 

t |a|<n, 

Z a ‘=td x j+d v j 

Now, note that the vector fields of the form td x j + d v j degenerate to d v j when evalu¬ 
ated att-0. Since from the Sobolev inequality, we have that 

||/(X,.)|| L oo (R n)< Y II 5 7(^--)IIz.1(R„), 

\a\<n 
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it follows that {T9) is striclty weaker than the inequality (10). In some sense, the fact 
that our method gives us a slightly worse estimate reflects its robustness and thus its 
appropriateness to deal with non-linear problems. 

Remark3.2. Estimates similar to fl8l holdfor the relativistic transport operator T m = 
\J 77? 2 + M 2 d t +£" =1 v‘d x i, where m — 0 for massless particle^ and m > 0 for massive 
particles. These estimates are presented in together with non-linear applications 
to the Vlasov-Nordstrom system. Interestingly, in the case of the relativistic transport 
operator, the estimates obtained have additional benefits. Decay estimates for rela¬ 
tivistic operators in the style of the Bardos-Degond estimate i 101 typically require the 
extra assumptions of compact support in v of the solution, while for the estimates 
obtained via the vector field method, only the finiteness ofthel} x v norms of the com¬ 
muted fields are required. See f4f. 

As is classical, derivatives enjoy better decay properties as follow^ 


Proposition 3.5 (Improved decay estimates for derivatives). For any sufficiently reg¬ 
ular function f defined on R" x R”, we have, for all reR" and all multi-index a 


lp(d“(/))|(x) < 


1 

(1 + t+ |x|)' i+|a| 


E ll^/lli 1 (R I xR„)i 

\fi\<,n+\a\,zP£.ff 


Similarly, for any sufficiently regular function f defined on R t xRjx R”, we have, for 
all t > 0 all xeU n and all teN, 


\p{d x d]{f))\{t,x) < 


(1 + f + |x|)" +|a| (1+ t) 


E \\z p f\\mu x > 

\f}\<n+T+\a\,Z^€.f^ 


Proof. The first part of the proposition is a consequence of Lemma [2T6l Lemma l2~4l 
and the global Klainerman Sobolev inequality 1 1 1 1 . 

The second part of the proposition follows similarly, using that 


tdtpif) — p{td t f) — p 


td t f-Y. vld v'f~ n f 


i =1 


with td t - Z” =1 V 'd v i being a linear combination of commuting vector fields as ex¬ 
plained in Remark lZ 1 1 □ 


Remark 3.3. In the above proposition, additional t derivatives yield only additional 
t decay and no improvement in terms of \x\ decay. This improvement can also be 
achieved assuming stronger decay in v of the initial data. More precisely, if f is a 
solution to T(f ) - 0 then 

dtf=-Y, v ' d x i f=Y J d x^- vl D- 

i =1 (=1 


Note now that iff is a solution, so is v 1 f, for any v l . Thus, for all i, f v d x i[-v l f)dv 
enjoys the additional decay (f + \x\) as stated in Proposition ^. 5\ provided that the L l 
norms of Z a {v' f) are finite for \a\ < n + 1. Iterating the procedure, we obtain that 
eachd t derivatives gives an additional decay of( t + x|). 

1 [n that case, the decay estimates are worse near the cone t = \x\, as for the wave equation, see again 

12 Note that however, quantities such as p[\d x f\) do not typically enjoy any additional decay. 
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For the applications to the Vlasov-Poisson system of this paper, we will also need 
an L 2 based Klainerman-Sobolev inequality to estimate VZ“</> pointwise, for <p a 
solution to the Poisson equation |2}. The estimate that we will used is contained in 
the following lemma. 


Lemma 3.2. For any sufficiently regular function ofyj defined on M”, we have, for all 


t> 0, 


ItgMI < 


1 

(1 + t + \x\) n ' 2 


E \\z a m\LHKY 

\a\<(n+2)l2,Z a ET s 


Proof. The proof is similar to that of Lemma [Q1 considering the function iy:y—> 
y/(x+ (t+ |x|)y) and replacing the L 1 Sobolev inequality with the l 2 Sobolev inequal¬ 
ity 


|i//(x)| 2 = |t/d0)| 2 < E f 

|a|<(n+2)/2*' B /i(°.l/2) 


\d a m 2 dy. 


□ 


4 Small data solutions of the Vlasov-Poisson system 

4.1 The norms 

The vector fields presented in Section [2~TI and the resulting decay estimates of the 
previous section will be sufficient to prove global existence of solutions to the Vlasov- 
Poisson system and derive their asymptotics for all dimension n > 4. On top of the 
boundedness of the L l norms of commuted fields Z“(/), which are needed in order 
to obtain pointwise decay from the vectorfield method, we will also need a little bit 
of additional integrability to prove the L p boundedness of the gradient of the com¬ 
muted potentials VZ a <p. 

With this in mind, for any n > 4, N e N and S > 0, let us consider, for any suffi¬ 
ciently regular function g of (x, v), the norm En,8 defined by 

£iv,<5[gl : = E H-Z (g r )lli 1 (R;xR£)+ E ll(l+M ) 2(1+5) Z (g)ll L i+«(RjxH;;)- 

|a|<IV,Z“ey l s “ l lalsiV.Z^ey^ 1 

As we shall see below, applying a similar strategy would fail in dimension 3 due 
to the lack of sufficiently strong decay. In order to close the estimates, it will thus be 
necessary to improve the commutation relation between our commutation vector 
fields and our perturbed transport operator. This led us to the introduction of mod¬ 
ified vector fields, denoted Y (and Y a for a combination of |a| such vector fields) 
below. Our main results are then similar to the n > 4 case, replacing the Z vector 
fields by the Y ones. In dimension 3, the norm En,8 will therefore be defined as 


£iV,5[g] 


E I! 

\a\<N,Y a £y$,s 


9 OIQ+3J 

+ E ll(l + M ) 20+ST Y (ghlyl+SfR^nja). 

|a|<iV,F“er'm! s 


( 20 ) 


The precise definitions of the modified vector fields Y and of the algebra J m ,s are 
given in Section [6721 
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Remark 4.1. In order to close our main estimates, we will not need to commute with 
any vector field containing t derivatives, i.e. commuting with vector fields in y s if 
n > 4 (respectively y m ,s if n = 3J will be sufficient. Commutations with vector fields 
containing t derivatives are of course usefull if one wants to obtain decay estimates 
of t-derivatives of p{f) and \7<p. In that case, one would simply modify the norms, 
replacing the algebra y s (respectively y m:S ) by the algebra y (respectively y m ). The 
interested reader can then verify that all the arguments below still hold. For these 
reasons, we will sometimes omit in the following section to specify whether the vector 
fields considered lie in y (respectively y m J or in y s (respectively y m , s )- 


4.2 The main results 


Our main results are the following 

Theorem 4.1. Let n > 3, 0 <8 < andN > 5n/2 + 2 ifn > 4, N > 14 ifn = 3. Then, 
there exists eq > 0 such that for all 0 < e < eq, if Fn,8 I/o I £ £, then the classical solution 
f{t,x,v ) o/CD-CD exists globally in time and satisfies the estimates 

1. Global bounds 

VreIR, E NiS [f{t)]<e. 

2. Space and time pointwise decay of averages of f 

£ 

for any multi-index a oforder\a\ < N— n, \p[Z a f){t,x)\ < 
as well as the improved decay estimates 
|p(d“/)(r, x) | ^ 


(l + m + |x|)' 


(l + |f| + |x|) n+l “ r 

3. Boundedness for L 1+s norms ofV 1 2 Z a tp 

for any multi-index a oforder\a\ <N, ||V 2 Z“^>|| L i + 5 (Rn) <e. 

4. Space and time decay of the potential and its derivatives 

for any multi-index a with \a\ < N- {3nl2 + 1), 

\Z a \I(p{t,x)\< ■ 

as well as the improved decay estimates 
\dfSI(p{t,x)\ < ■ 


f(n-2)/2 Q + | f| + |x|) 01/2 


t(n-2)l2 a + \ t \ + \ x \ ) nl2+\a\- 


Finally, all the constants in the above inequalities depend only on N,n,8. 

As explained above, when n > 4, the proof is easier and can be performed using 
only the commuting vector fields of the free transport operator while in the n = 3 
case, we will need to use modified vector fields. The modified vector fields approach 
would of course work also in the n > 4 cas(T~1but are not necessary there. In order to 
better explain the general framework, we will treat first the dimension n > 4 before 
turning to the proof in the n = 3 case. 

1 ’ A n interesting question left open by our work is to understand what happen in dimension 1 and 2 
where decay is even sparser than in dimension 3. We believe that at least in dimension 2, a carefull 

analysis using modified vector fields would lead to similar conclusions. We hope to treat this in future 
work. 
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5 Proof when n > 4 

In this section, we assume that n > 4, that the assumptions of Theorem l4.1l are sat¬ 
isfied for some initial data /o and we denote by / the classical solution of (T}-(3) 
arising from fo. 

5.1 Bootstrap assumption 

We will assume the following bootstrap assumption on the norm of the solution 
En, S [/( t)]. Let T > 0 be the largest time such that the following bounds hold 

Vfe [0, T], E N}S [f(t))<2e. (21) 

It follows from the smallness assumptions on En, sl/o] and a continuity argu¬ 
ment that T > 0. 


5.2 Immediate consequence of the bootstrap assumption 

Applying our decay estimate fl8t to \Z a f\ and \Z a f\ p (1 + v 2 )^ 12 as well as the im¬ 
proved decay estimates of Prouosition l3.5l we automatically obtain from our boot¬ 
strap assumption EIJ 


Lemma 5.1. For any multi-index a of order\a\ < N- n and for all t £ [0, T], 

Cn£ 


|p(|z a /|)(r,jc)| <----, 

|PU J 11 ' (l + |f| + |x|)" 


\p[\Z a f\ 1 +S a+v 2 )^){t,x) 


Cn£ 


1+5 


a+\t\+\ x \r 

for some constant C n > 0 depending on n, as well as the improved decay estimates, 

Cn£ 


( 22 ) 

(23) 


\p[d a f)(t,x)\< 


(l + |f| + |x|)” +l “ r 


5.2.1 Estimates on Z a ( <p ) 

From standard elliptic estimates, we can also bound an U’ norm of V 2 Z a <p 
Lemma 5.2. For any multi-index a with \a\ < N and for all t e [0, T\, 

||V Z 0(t)||^i+S(]jn) ^ Cjv,n,5£> 
where Cn,u,8 > 0 is a constant depending on N, n and S. 

Remark 5.1. There is obviously no difficulty in propagating higher IP norms ofV 2 Z a (p 
provided the initial data for f satisfy additional integrability and decay in v. We will 
not need them to close the estimates of our main theorem, which is why we did not 
assume the initial bounds on these U’ norms. 

Proof. Let p = 1 + <5. From the commuted equation for Z a <p and the Calderon- 
Zygmund inequality, we have 

HV 2 Z“(/)|| LP(R n ) < ||Z“(p(/))|| LP(H n ) , 

~ H IIp(|zV)|)lli P ( R «) 

|/3|<|a| 
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using Lemma [2741 Thus the bounds on V 2 Z a (f) follow if we can prove LP bounds on 
the p [\zPf\). For this, let us note that for any weight function j(id, we have, using 
the Holder inequality with IIp + 1/q - 1 

J \Z^(f)\dvj dx<J^J-^—jdv | dv^x^ v ) p ' q \Z^{f))\ p dvdx. 

Thus, we need to be integrable in v and we choose x(v) = (1 + \u\ 2 ) (S+n>12 . The 
lemma then follows from the bound on En, S [fit)] noting that with p = 1+5, we have 
piq- 8. 

□ 


Applying the Gagliardo-Niremberg inequality, we have immediately 
Corollary 5.1. Let q - For all te [0, T], 

||VZ“(/)(f)|| W ( R «) < C N} n }S e. 

Applying L 2 estimates for solutions to the Poisson equation (7) and the previous 
pointwise estimates on p (|Z“(/)|), we can also obtain L 2 decay estimates for VZ a (p 
provided |a| is not too large. 


Lemma 5.3. For all multi-index a such that\a\ < N-n 

\\VZ a <t>\\ L 2 {m < f( „_ 2)/2 - 

Proof. Multiply the Poisson equation satisfied by Z a <p by Z a (p and integrate by parts 
to obtain 


||VZ a (0)||2 2(RB) = - f Z a {(p)Z a {p{f))dx 

J X£lR n 


< [|Z“(0)|| 2n ||Z“(p(/))|[ 2n 

L7i=Z(R") r J L7nZ(R») 

Using the Gagliardo-Nirenberg inequality [ \y/ 1 1^ 2^ ^ < ||Vi//|| i 2( R n ) and Lemma l2.4l 
we obtain 


l|VZ a ((/))|| i 2 {R n ) < ||Z“(p(/))|| i _^ (Kn) 

< £ Mp(^(/))M i7 - (24) 

\p\<\a\ L (R } 

Since 

f dx 2 ^ f d[x,t) 2 < t 

Jxe R " (1 + | x | + f)n+2 ^ JCeR ” (| x / t | + 1)^+2 

the lemma now follows by estimating the right-hand side of J24) using the pointwise 
estimates m on p(Z^(/)). 


The previous lemma combined with the L 2 based Klainerman-Sobolev inequal¬ 
ity of Lemma l3.2l gives the pointwise estimates on V Z a (<p) claimed in the theorem. 

Corollary 5.2. For any multi-index \a\ < N—(?>nl2+ 1) and Z a e L 1 ^ 1 , wehaveforall 
fe [0, T], 

I I ^ 

l VZ ^1 ~ {\ + \x\ + t) nl2 f n -^ 12 ' 
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5.3 Improving the global bounds 

5.3.1 L 1 estimates of Z a f 

We first consider the L 1 estimate on Z“(/). From Lemma [23)1 we have for all multi¬ 
index \a\<N, and all t e [0, T], 

\\Z a f[t)\\ L i v < \\Z a f(0)\\ L i^+J^\\T ( i > (Z a (f))\\ L i ^ 

Thus, we only need to prove that the term below the integral is integrable in t. 

Now, from the commutation formula of Lemma l2.7l we know that 

l|70(Z“(/))|| i i <C N £ E K r ll|V i Z’ r 0V 1 ,Z^/|| L i 

|/3|<|a|-l,| r l+l/3|<|a| 

Note that we have so far no estimates on v derivatives of Z^f. To circumvent this 
difficulty, let us rewrite any v derivative as 

dyiZ^f) = (td x i +d v i)ZP{f) - td x ,z^f). 

Since both td x i + d v i and d x i are part of our algebra of commuting vector fields, we 
have 


\\T^Z a {fm L i v <C N {l + t) £ |C“ \\\V x zn(p)Z^f)\\ Llxv . (25) 

|/3|<|a|,|y|<|a| 

Irl+|j8|<|a|+l 


Now, since N > 5nl2 + 2 and [y| + \/3\ < \a\ + 1 < N + 1, we always have at least 
either |y| < IV- (3rc/2 + 1) or |/3| < N- n. 

Case 1: |y| < N- (3rc/2 + 1) 

In that case, we have access, thanks to Corollary [5j2] to the pointwise estimates, for 
all t e [0, T], 

\VzP(p{t,x)\<— Cn£ . . 
r (l + r )"" 1 

Thus, we get the estimates 


\w x zn(f))z^m L]cv 


< 


CjyC 

(1 + t) n ~ 2 


E ii^/iil 

l/3|<|a| 


l , 
x,v 


(26) 


and we see that if n > 4, then the error is integrable. 

Case 2: |/3| < N— n 

In that case, we estimate the error term as follows. First, 

|| V x ZH(f,)Z^f)\\ L i = [ [ \^ x Z^m\zhf)\dxdv 

x,v JxJv 

= J \V x Z r [(p)\p (|Z^(/)|) dx, 

since Z Y Up) is independent of v. Now applying the Holder inequality with 1 / p + 
IIq = 1, we obtain 

iWxzn^z^m^ < iiv,z^(0)n L j||p(iz^(/)i)|| 
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We would like to take q — 2, since the L 2 bounds are so easy to obtain for solutions to 
the Poisson equation. Using the pointwise bounds to estimate jp (|Z^(/)|) | , 2 , we 
would obtain 

||p(|Z'(/)l)|| 1|S c „«-“. 

On the other hand, we would get no decay a priori on ||V t Z r (0)11^2, since |y| is a 
priori too large to have access to decay estimates for the source term of the Poisson 
equation satisfied by Z r (0). With the extra weight of t in {25} , we see that if n = 4, 
we only get lit decay and we would get a logarithmic loss. 

To avoid this problem, we want to take q <2. Recalling the l/ 1 bounds of Corol- 
larv l5.ll we see that since <5 < , we have q - < 2 and thus, 


||p(|zV)|)|| i? < 


Cjye 

( 1 + t) nl2 ~ 1+,J ’ 


for some a > 0, which is now integrable in t. Putting everything together, we have 
obtain 


l|2>(Z“(/))|| z 


Cjve 


(1 + t) n ~ 2 


E n ,s [/(f)]- 


C N e 


C N e 2 


(1 + f) 


n/2-l+a 


d + f) 


n! 2-1+a ’ 


using the bootstrap assumptions to bound E n ,s [fit)]. 


5.3.2 v weighted I/’ estimates of Z a f 

Let p = 1 + 8 and q - 8{S + n). Recall from (8) the inequality, 

\\a+v 2 )« ,2 \z a fm p \\ LHR n xK) < ||a + ^ 2 )' ?/2 i^/(0)|P|| il(H „ xC) 

+ /j|(l + W 2 \Z a f\ v ~ X 2*(Z“/) | | il(R » xRS) ds 

+ J 0 II IZ a fjP(l + U 2 ) ( ‘ / - 1U2 d(P\\ Ll(R n xR n ) ds, 

< \\(i + v 2 ) ql2 \z a fm p \\ LHR>K) +h+i 2 , 

where 

a=£ lla+ v^iz^/r 1 T^z a f)\\ LHKxK) 

and 

h = /j||Z“/|P(l + , 2 )^ /2 d0|| il(R;>tR!;) . 

To estimate the error term I\, we proceed as in the previous section, replacing 
70 (Z“ /) using the commutation formula of Lemma [2771 and rewriting the terms of 
the form d l v Z& f as [td x i + d „) ZJHj) - td^zJ’ f. We are then left with error terms of 
the form 

(1 + f) f f |d x Z r 0||Z^/|(l + v 2 ) ql2 \Z a f\ p ~ l dxdv, (27) 

JxJv 

which need to be integrable in t. Using Young inequality, we have that 
f |Z^/|(1 + v 2 ) ql2 \Z a f\ p_1 dv < f \zPf\ p a + v 2 ) q/2 dv+ f |Z“/| P (1 + v 2 ) ql2 dv. 

J v J V J V 
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Note moreover that, as in the previous section, if \a\ < N - n, we have access to the 
pointwise estimates (23l . With this in mind, all the error terms coming from l\ can 
then be estimated as in the previous section. 

The estimates on the error term I 2 are easier and rely on the pointwise estimates 
of d(p of Lemma [5721 

5.4 Conclusions of the proof of the n > 4 case 

Thus, we have obtained 


En,S ( f ) - En,6 (0) + Cjv 


/' 


(1 + s) nl2+a 


ds, 


for some a > 0. Using the smallness assumption (0) < e, it follows that ( t) < 

e + Ce 2 < |e, provided e is sufficiently small, which improves our original bootstrap 
assumption ED and concludes the proof. 


6 The three dimensional case 

6.1 Strategy of the proof 

Repeating the previous argument in the n — 3 case would fail as seen from inequality 
(26}. Moreover, the non-linear terms do not seem to possess any special structure 
that would allow better decay such as the null condition for non-linear waves in 
3d, and it does not seem possible to close the estimates allowing the norms to grow 
slowly, using some hierarchy in the equations, as it sometimes happen for some sys¬ 
tem of non-linear evolution equations (for instance in 110)) - Thus, it seems that one 
is forced to try to improve the commutation relations so as to remove the most prob¬ 
lematic error terms. This will be done using modified vector fields. In hindsight, this 
strategy is reminiscent of the strategy of (2] where modified vector fields are con¬ 
structed by solving transport equations along null cones. 

To understand and motivate the definitions of the modified vector fields that we 
will use here, let us consider a solution / to the transport equation 'l]j, (/) = 0 and 
commute this equation with Z; = td x i + d v i. We obtain 


T,p(Zi(f)) = + fiV x (p.V v Zi(f) 

= Zdnf)) + piZi(y x <p.V v f)-^AZ i( p).V v f 

= ^(T 0 (/))~pf d xj Zi{<p).(td xj +d vJ )f + n £ d xj Zi{<p).td xi f 
7=1 7=1 

Thus, the error terms on the right-hand side are of two forms. The good terms are of 
the form d x Zi((p)Z' f , where the Z' are some of the commuting vector fields. These 
have enough decay so that they can be estimated as before. The bad terms are of the 
form td x Zj {(p)Z' f where, in view of the extra t factor, the previous arguments would 
lead to logarithmic growth. 

The aim of the modified vector fields will be to avoid the introduction of such 
bad terms. Note that on the other hand, commutations with vector fields such as d t 
or d x i would be better, because dd<p enjoys improved decay. As a consequence, we 
will only need to modify the homogeneous vector fields. 
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Let us consider, for all 1 < i < n, vector fields of the form 

n 

Yi = td x i +d yi - £ (. t,x , v)d xj , 

7=i 

where the coefficients ( t, x, v) are sufficiently regular functions to be specified be¬ 
low. 

Since d x j commute with the free transport operator, we have the commutation 
formula 


[T<p, YiKf) = - f i£,d x ,Zrt<l))Z j (f )+iiY J ®[d xi {V x <p).V v f (28) 

7=1 7=1 

- t T ( p{<&l)d x j f + p £ d xj Zi«p) td xJ f 
7=1 7=1 

where Zj = td x j + d v j. The first term on the right-hand can be handled as before as 
it does not have v derivatives leading to the extra power of t. Note moreover that 
since td x j is part of the algebra of macroscopic commuting vector field, the second 
term can be rewritten as 

Ztft- l \ x Zj<p\ v f 
7=i 

and thus is expected to be integrable by the previous arguments provided are 
uniformly bounded by any power of t strictly less 1. 

The key idea is then to choose appropriately the 0 J . so as to be able to cancel the 

last two terms in (28). For this, we impose that each is obtained as the unique 
solution to the inhomogeneous transport equation 


70(<J>j) = iitd xi Zi{(p ) 


with 0 initial data. Note that the right-hand side of this equation only decay like 1 It, 
so we expect <T j to actually grow logarithmically in t. 

Assuming that this holds, we see that all the error terms in l28l are now inte¬ 
grable and thus, we should expect to control our norms EN,slf{t)] provided all the 
original commutation vector fields Z are replaced by their modified form Y, hence 
the definition of the norm in the 3d case given by l l20l . A few difficulties remain 

1. Since the (and in fact all the coefficients involved in the contruction of the 
modified vector fields) are growing logarithmically, it is apriori not clear how 
to exploit the new energies obtained after commutation with modified vector 
fields to obtain pointwise decay. The idea is again to rewrite the extra terms 
such as fJ)^ d x j (/) in the form 


if) - -f {td x j + d v j ) (/) - % V J (/) 

and to use an integration by parts in v for the last term to push the v deriva¬ 
tives on the ([>1 coefficient. See Section lOl below. 
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2. The Poisson equation A <p = p{f) cannot be commuted with modified vector 
fields, since the coefficients in the modified vector fields depend on u, while <p 
and p(/) are macroscopic quantities depending only on ( Thus, we keep 
commuting this equation with non-modified vector fields. Quantities such as 
Z a p{f) are then rewritten as p(Y a {f )) plus error terms. The structure of these 
error terms is the subject of Lemma l6.3l 

3. As is seen from the statement of Lemma [6731 some of the error terms will de¬ 
pend on Y a {(p), where <p is a coefficient obtained by solving a transport equa¬ 
tion of the form 2^ (cp) - td x Z(p and Y a is a composition of | a | modified vector 
fields. Commuting the transport equation satisfied by cp by [a| vector fields, 
we see that we need to control tY a d x Z{(p). This poses a problem at the top 
order, when \a\ - N, since it looks as if one needs to control dZ a Z[(p), which 
would a priori require to commute the Poisson equation by N +1 vector fields, 
and therefore would forbid us to close the estimates. On the other hand, us¬ 
ing directly the Poisson equation satisfied by Z“ (</>), we may hope to control 
d 2 Z“ (</>). To exploit this fact and close the top order estimates, we devote some 
time to describe the structure of the top order terms of Lemma l6.3l in Lemma 
E2] 

4. Finally, numerous terms of the form Y a {<p)Y^{f), where <p is a coefficient ob¬ 
tained by solving a transport equation as above, will appear in the equations. 
This creates another difficulty to close the top order estimates, since we do not 
have access to pointwise estimates on Y a ((p) when |a| is large, nor do we have 
access to L p x v estimates on Y a ((p) because the source terms in their transport 
equations are not integrable in v. We will circumvent this difficulty by consid¬ 
ering directly a transport equation satisfied by the product Y a {<p)Y^{f). See 
Section l6.5| below. 

We now turn to the details of the proof. 

6.2 Definitions and properties of the modified vector fields algebra 

As explained above, our new algebra of microscopic vector fields will consist in 

1. Standard translations d t , d x i, 

2. Modified uniform motions Yj := td x i + d v i — £? =1 ^ k B x k, 

3. Modified rotations Of . + QV , wf .d r k, 

4. Modified scaling in space S x + S v -X£ =1 ° k d x k, 

5. Modified scaling in space and time td t + Z" =1 x l d x i — £? =1 B k d x k, 

where the coefficients <£> k , u> k ., a k , 6 k are solutions of the following inhomogeneous 


equations with 0 data at t — 0 


Z*(of) = ptd x k[Zi{ip)\, 
T,f,(ciJ k j) = ptd x k 
r 0 (CT fc ) = ptd x k[S x (<p)-2(p\, 



T<p(O k ) = ptd x k td t + Y J x l d x t (</>) 


2 — 1 
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6.2.1 Further notations 


Since the sign of /i will play no role in the analysis to come, we will asume without 
loss of generality that p = 1 in the rest of this article, to symplify the notation. 

We will denote by Jl the set of all the coefficients a) k o k , 6 k and by <p a 
generic coefficient among them. Similarly to the set y and y s , we define the sets y m 
and jm iS where y m is the set of all modified vector fields (including the translations) 
and Y m ,x is the set of all modified vector fields minus the time translation and the 
space-time scaling. 

If Z is an original, non-modihed vector held, we will sometimes write schemat¬ 
ically Y — Z + (pd x to denote the associated modified vector held, where by conven- 

n 

tion <pd x - 0 is Z is any of the translations, cpd x = E if Z is one of the uniform 

fc=t 

motions and similarly for the other vector helds. 

Finally, we will use the notation P((p) to denote a function depending on all the 
coefficients in y m or y m , s . 


6.2.2 Improved commutation formulae 

With these dehnitions, we now have the following improved commutation formulJ^l 

Lemma 6.1. For any Y e y m and any sufficiently regular function g of [t , x, v), [, Y](g) 
can be written as a linear combination with constant coefficients of terms of the form 
d x Z(0)F(g), <pd x Z((p)Y{g) or <p 2 d x Z((p)Y{g) where ip e M, <p 2 denotes a generic 
product of two coefficients in dl, ZeT and Y e y m . 

Proof. If Y is a translation, for instance Y = d x k, then 


[T<p,Y]{g) 


-{d x kV x (t)).V v g, 


~T. d x‘ d x k( P- td x ig + d v ig- Z,^id x jg 

i=t i =i 


E d x< d xkf>- td x t g - E ® • d x i S 
i'=t ]=i 


-'L^xidxktp.Yiig) 


i= 1 


E d x‘ tdxt(p.d x i g - E E ® • d x id x k (0) ,d x j g, 


i=l 


i=1 j =1 


which is of the desired form since td x k £ T. If Y = Y[ = td x i + d v , - Z" = i ^xi > then 

it follows from l28l and the dehnition of the coefficients that 


\T, P , Yffig) = - E djZtWZjlg) + E ®^x/(V*0).V„g. (29) 

7=1 7=1 


where Z; (0) = td x i (p and Zj (g) = td x ig + d vl g. Since Z ; -(g) = Yjig) +Z'Li ® l jd x iig), 
the hrst term term on the right-hand side of J29) is of the desired form while for the 

14 Once again, in all the formulae of this section, we assume that <l> is a sufficiently regular function of 
(t, x) defined on [0,7' x R 3 and we will eventually control the regularity of </> through a bootstrap argu- 
ment. 
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last one, we have 


7 = 1 


= Y J Y,®\ d xi d x* , P- d vk8 

j=\k=l 

= tt^i d xid x ^.Y k {g) 
j-lk-l 

n n n ' 

-E E®/ a */ fl **0- - 

j=ifc=i ;=i 

= E E ^ J i d xJ d x k( P- Y ic8 

j= ifc=i 

-EL + E E it^Mdxitfx^xig’ 

j=lk=l ]=lk=ll=l 


which is of the desired form. The commutation with the other modified vector fields 
can be computed similarly. Note that in the case of the spatial vector field, we have 
[T^,, S x + S 1 ’] = -SJ x S x ((f)).Vv + 2 V x (p.\7 v , which explains the extra term in the defini¬ 
tion of the coefficients cr k compared to the other vector fields. □ 

If we apply several times the previous formula in order to compute [r 0 , Y a ], 
many products of the form Y p (qf)Y v {<p') will appear. In order to simplify the pre¬ 
sentation below, it will be usefull to use the following definition. 

Definition 6.1. We will say thatP{(p) is a multilinear form of degree d and signature 
less than k ifP{(p ) is of the form 

pm= E c p n yPi W’ 

pel d , j=h-,n, 
p=Cpi,--,Pd) 


where I denotes the set of all multi-indices (thus pj is a multi-index for each j), for 

d 

each p in the above formula E IP/'I - fc and where the C p are constants. 

i=i 

From Lemma [PI we now obtain 
Lemma 6.2. For any multi-index a, we have 

[r 0 ,7“] = '“£ 1 f £ P“jj<p)d xi zn(l>)YP, (30) 

d= 0 1=1 |y|<|a|, 

\P\<\a\ 


where the (ip) are multilinear forms of degree d and signature less than k such 
thatk< \a\ - 1 andk+ |y| + \ f\ < |a| + 1. 

Proof. This is a classical proof by induction for which we will just sketch the details. 
Lemma lim shows that l.’-l()l holds when \a\ = 1. Assume that l30i holds for some 
multi-index a and let Y e y be an arbitrary modified vector field. We have 

[2>, YY a ] = [r 0 , Y] Y a + F[T 0 , Y a ]. 
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One easily see that the first term on the right-hand side has the correct form using 
Lemma I6T1 The second term will generate three types of terms. The terms of the 
form 

Y ( pa 4p] d x‘i zr w) Yp 

are of the correct form, the multilinear form being of the same degree and its signa¬ 
ture being increased by 1 at most. For the terms of the form 

we recall that Y is schematically of the form Z + (pd x . Thus, 

Pd'yB Y (d x t Z r {(p )) Y@ = £ P'd^Z^mY?’, 

ir'islrl+i 

where P', are multilinear forms of degree at most d +1 and we have not changed the 
signature, so that these terms are of the desired form. 

Finally, the last terms are of the form P l *^,d x i (ZJ (0)) YY@, which clearly satis¬ 
fied the required properties. □ 

As explained above, we also need to revisit our commutation relations for the 
Poisson equation satisfied by the potential <p. Contrary to /, we cannot commute 
with a modified vector field, because the coefficients of the modified vector fields 
depend on v, while (p is a macroscopic quantity and depends only on ( t, x). Thus, 
we keep commuting with Z“. We have 

Lemma 6.3. Let g be a sufficiently regular function of[t,x, v) and let(p g solves 


A</> g = p(g). 

For any multi-index a with \a\ < N, we have, for all 0 < t < T, 
\a\ |o|+l i , , l«l 


p(Z“(g))=EE X -jP{P a d iWY fi {g))+t. E p[Qd^ x( p)Y^g)) (31) 

1-1 \R\<r\n,\ V H \ R\<r\rv\ V J 


and 


y=l d =1 ipi<|ct| 


l«l |a|+l 1 , . „ . \a\ , 

AZ“(0)=EE E -p(^(<p)^(g))+E E p[QdpMY P lg)}, (32) 


d=0 |/3|<|a| 


|o| 


7=1 d= 1 |/3|<|a| tJ 


d= 0|pi<|a| 


where P{cp) and P^j {cp) are multilinear forms of degree d and signature less than 
k satisfying 

Ic < \a\, k+\fi\< \a\ 

and where the Q^{d x cp) and Q^{d x cp) are all multilinear forms of degree d of the 
form 

Q{d x (p) = E C P FI Y Pi ^ x rjcp), (33) 

pe/ d , j=h~,n, 

p=(pi,..,p d ) 

d 

where the C p are constants, 1 < rj < 3, and such that k' := E P / 1 satisfies 

f=i 

k' <|a|-l, d + k' + \fi\ < \a\. 
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Proof. First, recall that if Z £ F, then 

A Z{(pg) = Z\<p g + d.zk.(pg, 

where dz — 0 unless Z is one of the two scaling vector fields, in which case d/ - 2, 
and 


Z(p(g)) = p(Z(g)) + c z p{g), 

where cz — 0 unless Z is the spatial scaling vector field, in which case cz — 3. Since 
A <p g — p (g) , it follows that l3D implies I .'-> 2 1 . 

In the next lines of computations, given Zey and Y the modified vector field 
corresponding to Z, we will use the schematic notations Y — Z — cpd x and Z - Y + 
ipd x instead of any of the lengthy formulae given at the beginning of Section 16.21 
such as Yi — td x i + d v i — L£ = 1 0^d x k. We will also use the notation td x + d v - <pd x to 
denote a generic vector field among the Y t , the letter Y' to denote a generic modified 
vector field and the letter <p' to denote a generic coefficient belonging to .M. 

We now compute, for any Z £ y 



/ {Z + (pd x -(pd x )[g)dv 
J V 

/ Y{g)dv- / ipdxgdv 

Jv Jv 

/ Y{g)dv- ^(td x g + d v g-(pd x g-d v g + cpd x g)dv 
JV J V t 

f Y{g)dv-~ f (p[Y\g)+cpd x g)dv+ f ^ d v gdv 
Jv t Jv Jvt 


The first and second terms on the right-hand side of the last line have the correct 
forms. For the last term, we integrate by parts in v 



- f d v cpgdv 
t Jv 

7 f [td x + d v -(p'd x -td x + (p'd x ){cp)gdv 
t Jv 

7 f [Y 1 + <p'd x ) [<p)gdv- f d x (cp)gdv, 
t J v J V 


where now all terms are of the correct forms. This prove J3T) when |a| = 1. We 
now assume that J3T} is true for some a. Let Z be a non-modified vector field. Us¬ 
ing again that Zp(Z“(g)) = p(ZZ“(g)) + Czf>(Z a (g)), we only need to prove that 

Zp(Z“(g)) is of the correct form. Assume first that Z contains no f-derivativtf^l, 

n 

i.e. Z i=- d t and Z ± td t + ^ x l d x i- Using the induction hypothesis and writing Y = 

i =1 

Z + cpd x to denote the associated modified vector field, we have 

15 Recall also that commuting with the vector fields containing f-derivatives is not necessary to prove 
Theorem l4. 1 l and is only usefull if one wants to obtain improved decay of t- derivatives. 
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Zp(Z“(g)) = z 


' |a| |a|+l 1 , . , \a\ 

LEE (g) +E E p 

Vj=l rf=l iPls|a| 4 J 1 


EE E ±jp(z[p%mY p ig)])+'E E p(z[Q&dx<p)Y p w]) 

]=1 d =1 |/3|<|a| J ’ 


d=HPI<|a| 
|a| 


d=HPIs|a| 


+ c z p{Z a {g)). 


(34) 


If now Z contains t- derivatives, we would get extra terms in 1341 which arise when 
d t hits the -y factors. Since these extra terms are all of the correct forms, so we only 
need to analyse the terms on the right-hand side of J34). 

The last term in (341 has already the right form. Next, replacing Z by Y + <pd x in 
the terms p ( Z [ (<p) Y& (g) ] j, one easily see that they also have the desired form. 

For the terms p \z [ Q“^ (d x <p) Y‘ 0 (g) ] ], we have 

p[z[Q^d x (pW^g)]) = p({Y + cpd x )\Q^d x cf»Y^g)]) 

= p(Y[Q a dfi ( d x cp ) Y? (g) ]) + p (cpd x [ Q a dfj {d x (f» YP (g) ]) 

The first term on the right-hand side is easily seen to have the correct form. For the 
second term, we write cpd x - j ( td x +d v - q>’d x + q>’d x - d„) = jY’ - jd v + -j-d x so 
that 

p[<P d x {Qdp(d X (p)Y ft {g)^j = ~ t p[{(pY' + <p<p'd x ) \Q dp (d x <p)Y p (g)^ 

+ -p [d v <p\Q%p(d x <p)Y p (g)] ] 

using an integration by parts in v. The first term on the right-hand side has now the 
right-form. For the second term, we again write d v cp = ( td x + d„ - (p'd x ) <p - td x <p + 
q)'d x <p so that 

- t p[d»<p\Qdpld x <P)Y p lg)]} = - t p[Y'l<p)[QSpld x fp)Y P lg)]} 

-p[d x l<P)[Q^d x v)Y p (g)]) 

+ ~ t p[(p'd x ((p) [Q d pid x ip)Y p (g)^ 

where all terms now have the correct form. □ 


Recall that we do not hope to have any good estimate on F“(<p) if |a| = N, since 
its transport equation would then contain a source term of the form VZ“Z</> and we 
only hope to have estimates for VZ^0 up to fi — N. On the other hand, using directly 
the Poisson equation satisfied bv Z a (d>), we will have good estimates on d x d x Z a Uf>). 
To take advantage of that fact 1 ! we will need later the following technical lemma, 
which improves upon Lemma [631 bv describing the structure of the P multilinear 
forms a little further. 

‘■Note that these difficulies arise only at the top order |a| = N. An alternative to the approach taken 
here would be to allow for the top order estimates to grow slightly in t. 
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Lemma 6.4. With the notations of Lemma \6.3\ the multilinear forms Up) and 
P^’p ((p) can be written as 


P di [(t>) = P ' [<P) + E E E &i< P «P)d x ‘Y Tli M'{(p), 


q)Z.M i =1 


/=1 


= Pi(^)+ E E E* 3 V0> a *‘ yW w. 


i=l 




where 


1. P\ and P\ are multilinear forms of degree d and signature less than k satisfying 
k<\a\-l,k+ \f \ < \a\, 


2. Yi are the modified uniform motions defined at the beginning ofSection \6.2\ i.e. 

Yi = td x i+d v i- Z n k=l ® k id x k, 

3- \pi,<p\ = |a|-l, Pt,cp = \a\- 1, \rn, Vl(p ’\ < |a|-l, I f] iiipi(p f\<\a\-l, 

4. Ci iip , C (> are constants and & > i (p {cp'), SPiqyfip) are polynomial of degree at most 
\a\ in thetp’ eM. 


Proof. This is an easy proof by induction. The case \ a\ — 1 has already been proven 
in the proof of the previous lemma. Assume that the statement of this lemma holds 
for some a and let Z e T. As before, it is sufficient to consider only Z (p(Z“)). We 
will only be interested in the top order terms, that is to say terms containing Y' 1 Up) 
with rj = | a | + 1. Note that they can only be generated by applying a vector field to 
terms containing Y r> Up) with rj' - \a\. 

Using the induction hypothesis, the top order terms coming from the P multi¬ 
linear forms will give terms of the form 


E E C it(fl Z[YiY Pi -’n<p)\ (35) 

cpeM i= 1 


and 

E 06 ) 

ip.cp'eM i=l 

For the first type of terms, we have 

Z[Y*WM<p)] = lY + <p"d x ) [YiYP^Up)] = YYiYP^Up)+ip"d x YiYPi-ncp). 

Now the second term on the right-hand side has the right structure. The first term 
also has the right structure, since [Y, Yf can be written as a linear combinations of 
terms of the form d x , the modified uniform motions Yj and <p"d x and Y'Up)d x . 

The terms of the type l36l can be treated similarly, the only dangerous terms 
being of the form cp'Yd x i Y^'w' (<p), where the Y and d x t can be commuted up to 
lower order terms. 

The top order terms coming from the Q multi-linear forms will give terms of the 
form 

p(z[r»(d* V )F^(g)]) = p[{Y + cp"d x )\YP{d x cp)Yf i {g)\) 

= p (F [ Y* [d x cp] Y?{g) ]) + p [(p"d x [ YP [d x cp) Y? (g) ] ) , 
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where \rj\ = \a\ -1. Now the first term on the right-hand side will contribute only the 
Q forms so it can be ignored here. While the second term, repeating the argument 
of the previous lemma, gives the following contribution to the P forms 

- t p(ncp)"Y , Hd x (p)Yf i (g)) 

which, since \tj\ < \a\ - 1 is not of top order and 

as well as 

~ t p(cp'cp"d x [Y\d x( p)Yf i {g)\ ) j = - t p[cp , cp"d x [Y^d x cp)]Yf , {g) ) j 

+ ip(<pV'F I '(d^)d i -[7^(g)]) (37) 

which are all of top order. Again, we use that d x essentially commutes with any 
modified vector field up to lower order terms, to put these last terms in the right 
form. 

□ 


The following commutation property will also be usefull later. 
Lemma 6.5. For any multi-index a and any 1 < i < 3, we have 

\d x i,Y a \=f, E p -iS d x <P)Y P d x j, 

d =0 |j6|<|a|-l, ’ 

1<;<3 


where the id x cp) are multilinear forms of degree d of the form l33l with a signa¬ 

ture less than k such that k+\f \ <\a\ - i . 

Proof. We prove the |a| = 1 case, the general case following by an easy induction 
argument. Let Y be modified vector held. We will write schematically Y — Z + <pd x , 
where Z is a non-modihed vector held and cpd x stands for a linear combination of 
products of some (p e Jl and some d x k. For simplicity, assume that Z commutes 
with d x i (the cases where Z do not commute with d x , can be treated similarly, since 
the resulting error terms do not involve any terms depending on the coefficients in 
,M.) An easy computation then shows that 


[d x i,Y] = [d x i,(pd x ] = d x i{(p)d x , 

which is of the desired form. □ 

Finally, let us also remark that 
Lemma 6.6. For any multi-index a, 

1 |a| 

Y a V<p — Z a \/(f>+ - ]T P%(iip) z ^<p, 

1 d =1 p 

where (<p) are multilinear forms of degree d and signature less than k such that 
k < |a| - 1 and k+ \f \ < \a\. 
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Proof. We only do the |a| = 1 case, the rest of the proof being an easy induction. Let 
Y' be a modified vector field and write schematically Y' — Z' + (pd x . We have 

7'V0 = (Z r + cpd x ) V</> = Z'V<p + ^ [ td x ]V(p, 

which is of the correct form since td x ef. □ 


6.3 Bootstrap assumptions and beginning of the proof 

Let now / be a solution to the Vlasov-Poisson system in dimension n — 3 such that 
the hypotheses of Theorem l4. 1 l are satisified. 

We consider the following bootstrap assumptions. Let T > 0 be the largest time 
so that, for all t e [0, T] and all x e l 3 , 


1 . 


E Ni5 {f{ t)]<2e, 


(38) 


2 . 


For all 0 < S' < 8, there exists a Cy > 0 such that for all multi-index a with 
\a\ < N, 

\\V 2 Z a <pm\ L i +S '<C S 'e 112 . (39) 


3. For all multi-index a with \a\ < IV- (9/2 + 1), 

e 1/2 

|VZ“</>(f,x)|<--- (40) 

1 +t 

and 

£ II2 

\Y a V(p{t,x)\ < -- (41) 

1 +1 

4. For all multi-index a with |a| < N- (9/2 + 2), we have for all cp e ,M, 

| Y a (cp) ( t, x, v) | < e 1 ' 2 11 + log(l + t) |, (42) 


5. For all (p e Jl, all 1 < i < 3 and all multi-index \a\ < N- (9/2 + 3), 

\d xi Y a {q)){t,x,v) \ <e m , (43) 


It follows from the initial data assumption and standard arguments that T > 0. In 
the rest of the proof, we will try to improve each of the above assumptions, which 
would show that T = +oo. 

Note that in view of Lemma [63] j43) is equivalent to 
\Y a d x t{<p)V,x, v)\<e 1 ' 2 

and we will switch freely between the two in the rest of the article. Note moreover 
than in view of the Gagliardo- Niremberg inequality, assumption f39t immediately 
implies that 

||VZ a 0||i9 < Cq£ 112 , (44) 

with q = 3/2 < ^-s' ^ < in view of the definition of S. In particular, q can be taken 
as close to 3/2 as wanted. 

Finally, the notation A< B used in 1401 . II 1 1 . l42i . (43) stand for A < C^^.nB 
where Cjv,5,n is a constant depending only on N,n,8 (with here n - 3). We will 
eventually improve each of these inequalities by replacing the e 112 on the right-hand 
sides of l40l , 07}, (42) , (43) by e, choosing e sufficiently small to absorb the CN,s, n ■ 
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6.4 Klainerman-Sobolev inequalities with modified vector fields 

Using the bootstrap assumptions above, we have 
Proposition 6.1. For all multi-index\a\ < N — 3, 


p{\Y a nt)\) < 

< 


i 

(l + f+|x |) 3 


E 11^/coiiii^xRi, 

|j8|<|a|+3 


1 

(l+f+|jc|) 3 


E II^/(0IIl1( R 3 xK 3 ) . 

\P\ <N 


Proof. Similarly to the proof of Proposition [373] let us fix (t,x) eRx R 3 and let yr be 
defined by yr : (0,1/2) 3 y — p (| F“/|) [t, x+ [t + \x\)y). We fixd' = ^ and apply a 

Id Sobolev inequality 

p{\Y a f\)[t,x) < f (\d yi yr\+\y/\){y lt 0,0)dyi, 

J\yi\<S’ 1 ' 2 


where as before 


d yi y/[y) = {t+\x\)d Xl p{\Y a f\)[t,x+(t+\x\)y) 

= t f d Xl (\Y a f\)(t,x+{t+\x\)y, v)dv + \x\ f d Xl [\Y a f\ (t,x+ {t + \x\)y, i/)) dv. 

J V J V 


Now, 


tf d Xl [\Y a f\)dv = f 
Jv Jv 


td Xl +d Vl - £ <b{d x ] + £ <t>{d xJ 
«=1 7=1 


(|F“/|) 


dv 


= f Yi[\Y a f\)dv+ f £®[d xi {\Y a f\)dv. 

Jv Jvj~l 


The first term on the right-hand side is simply estimated by 

If Y\ (|F“/|) dfl < f \Y\Y a f\dv. 

I Jv I Jv 

For the second term, we again try to force the apparition of our modified vector 
fields 


f <p[d xi (I Y a f\) dv = f ^ td x j + d v j - £ <s,)d x . 

Jv Jv t \ k=1 


(I Y a f\)dv 


-f v Y[dvJ-Zfj d x f\Y a f\)dv 
f ^ Yj{\Y a f\)dv 

Jv I 

-Xy(^-E®^) [\Y a f\)dv. 


The first term on the right-hand side can then be estimated as above, using that 

_ 

| -j -1 is uniformly bounded from the bootstrap assumptions 1 1421 . For the remainder 
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terms, we first note than in view of the bootstrap assumptions (42) , the terms of the 
form 

can be estimated by 

f \dA\Y a f\)\dv. 

Jv 

For the last type of terms, we integrate by parts in u 


r (Er r ■ (!)■' 

-±d vJ (lY“fl)dv = - ——— [\Y a f\] d 

Jv t Jv t 


V. 


We now rewrite d !jJ as 


dv>®[ = 


‘ td xi + d vi - £ ofaJ - td xJ - £ ^d xk | ^ 
I k= 1 


k= 1 


= Yj{ ®i)- 


td xJ-£^d xk U{. 


k=l 


The first term only grow like 11 + log(l + t) | according to (42) and this growth can be 
absorbed thanks to the 1/ f factor. For the second term, using (43) and (42) 


I td x j (0{)| + \^d xk (0{) I < e m t + e 1/2 | 1 + log(l + r) I 


and again we can absorb the growth using the 1 /1 factor. 
Putting everything together we have obtained that 


pi\Y a {f)\)(t,x) < f [ [\YY a (fl\ + \Y a [fl\)[t,x+{t + \x\){yi,0,(»)dvdy 1 . 

J\yi\<S ' 112 Jv 

The remaining of the proof follows as in the proof of 13.3) , repeating the previous 
arguments for each of the variables and applying the usual change of coordinates. 

□ 


6.5 Estimates on products of type Y a [tp)Y^ (/) 

Due to the form of the commutators of Section l6.2.2l we will need to estimate terms 
of the form Y a {cp)Y^ (/). When a is sufficiently small, we will have access to point- 
wise estimates on Y a (<p) so there is no difficulty. When a is large, say a = N, we 
have for the moment no estimate on Y a {ip) and we can certainly not hope to prove 
pointwise estimates for these quantities, because, in view of the transport equations 
satisfied by the coefficients tp, these estimates would in turn require pointwise esti¬ 
mates on VF“ ((/>), and these estimates do not hold at the top order. Instead, we will 
prove directly estimates on the products Y a {ip)Y^(f), taking advantage of the fact 
that YP{f) are integrable in v. More precisely, 

Proposition 6.2. Let cr > 0. Then, there exists aC a > 0 such that for any multi-indices 
a,f, | cr | < IV— 1, |/3| < 9/2+3 and any <p e M, we have for all t e [0, T] and all 1 < i < 3, 
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\\Y a {cp)WY^f)m L H« x m v ) ^ C a a + t) a e, 
\\YiY a ((p){t)Y^f){t)\\ L i (KxXUv) < Cail + tfe, 
\\d xi Y a mt)Y^f)m L HR x m v ) s Cad + fl^e, 

Proof. Let <7 > 0 and let a be a multi-index satisfying |a| < ;V and such that if \a\ - N 
then Y a - Yj Y a ' or Y a - d x j Y a ' with \ a'\ = N- 1. 

We have 

T^[Y a [<p)Y p [f)] = r 0 (y a (v»)) f^(/) + y“(^)r 0 (y^(/)) = 7, +/ 2 , 

where 7i = T<p[Y a {ip )) F^(/) and 7 2 = Y a {/)). In view of Lemma 12.91 it 

suffices to show that 

l|7l,7 2 || L l ^ < Ccrfl + t) l ~ a . 

1. Estimates on I\ 

We have 


h = [^,F a K V )F^(/) + F a [r 0 ( V )]F^(/) 

= A, I + fl,2, 

with 7 U = {T^Y^WYHf) and 7 i , 2 = Y a [T^<p)]YHf). For 7 U , we use the 
commutation formula t30l 

d=0 i= 1 |yl<|a|, 

|j)|<|a| 


where satisfies the requirement of Lemma Rx2l in particular, it has signa¬ 
ture less than k such that k < \a\ - 1 and k+ \j\ + I 77 I < |a| + 1 . 


Case 1: |y| < IV-(9/2 + 1). 

It then follows from the bootstrap assumption l40l . that we have the pointwise 
estimate 

£ 1I2 

K<z r w\<— . 

Moreover, since k + \r/\ < \a\ + 1 < N + 1, either k < N- (9/2 + 2), and we have 
access to the pointwise estimate 

|p“^(<p)|<c d/2 (l + log(l + r)) d 

or k > N- (9/2 + 2). In this case, we have \ri\<N+l-k<N - (9/2 + 2), since 
N> 14, so that we now have access to pointwise estimates on Y'Hcp). Note also 
that since k < N + 1, there is at most one factor in each of products of Y p i Up) 
in the decomposition of W) for which we do not have access to pointwise 
estimates. In conclusion, it follows that we have an estimate of the form 


Pjj v m X iZi'{ ( p)Y\c P )\ 


1/2 (1 +log(l + D) rf 

1 + t 2 


E \Y p m- 

|p|<|a|-l 
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Case 2: y > N— (9/2+ 1). 


Then, fc+|? 7 l < IV+l-|y| < N- (9/2+2) since N > 14 and we can bound (<p) 

and Y^Up) pointwise. Since |/J < 9/2 + 3, we have access to pointwise bound 
on p(|Z^(/)|), so that we can estimate 

| |P“^ md x i Zr (0) Y* ( V ) (/) 11 4 < 

e 1/2 (1 + log(l + t)) d+1 \\d x i I l i? I Ip (l Y? (/) |) | | L j, 

where IIq +11 p - 1. Taking c/ as in J 11 1 with S' as small as needed (depending 
only on er), we obtain, using the pointwise estimates of Proposition l 6 .1 1 that 

llp(irV>l)ll I 5SC„_i_ 

so that 

^ e 2 (1 ■ 

which, assuming a < 1 , is integrable in t. 

We now turn to the estimates on I\ : 2 - Recall that we have T^Up) - td x i Z (</>) for 
some Z and some x 1 , unless Z is the spatial scaling vector held, in which case 
T(p{cp) = td x i (Z{(p) - 2(f)). Since the extra term can be handled similarly, we 
will only treat the case of the non spatial scaling vector fields below. Applying 
Lemma [6761 we have 

Y a { T^cp)) = tY a d x iZ{<p) 

1 |a| 

= tZ a d x ,Z(<l» + t-Y P a d A<p)Z*d x iZW), 

1 d =1 
\a\ 

= tZ a d x iZ[(f>)+Y J Pdr ] mZ rl d x iZ{(p), (45) 

d =1 

where Up) is a multi-linear form of degree d and signature less than k with 

k< |a| - 1 and k+ \q\ < \a\. 

Assume Hrst that |a| < N- 1. For the first term on the right-hand side of (45), 
we have 

\tz a a x ,zm<t y \d X iZ^zm. 

\rt\<\a\ 

Since \a\ < N—l, we have I 77 I + I < N in the above sum. Thus, \tZ a d x iZ{(p)\\Z^[f)\ 
can be estimated as before using the Holder inequality, pointwise estimates 
on p (| (/) |) and the estimate 1 +TI . 

If now |a| = N, then by assumption, Y a - Y^ Y a ' or Y a - d x j Y a ' with \a'\ = 
N- 1, so that Z a = td x iZ a ' or Z a = d x iZ a '. Thus, 


\tZ a d x iZ{(p)\ < 
< 


ta+t) Y 1 dxdjz'zm, 

fd+f) Y \&x zV w\- 
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We can now estimate \tZ a d x iZ((p)\p(\Z^{f)\) using Holder inequality with 
p = = 3 - = 1 + o' with a 1 = 7 ^ > 0 assuming a < 3 and q - j- ; . Recall that 
Wd^-Z^ (<j))\\ L p is bounded thanks to the boostrap assumption 1391 provided a 
is sufficiently small. Moreover, using the pointwise estimate on p(|Z^(/)|), we 
have 


llp(l^(/)|)||^ (R 3 ) <C(l+t) 3 - ,T . 

Whether \a\ < 1 or \a\ = N, the above estimates gives 

r e 3 ' 2 

\tZ a VZm\Z f) [f)\dxdv<C <7 - --y— 

J x,v U + t) 

which, after integration, gives rise to the t a growth in the statement of the 
proposition. 

For the second term on the right-hand side of J45), we have either |p| + 1 < 
N- (9/2 + 1), in which case, we have access to the pointwise estimate 

e 112 

ia -' z ’ zwls iITtF 

and thus, 

p a d md X iZ^zmz^f)\< e ( \V!l 1+r)) E \Y p m\zhf)\, 

where we have used the fact that there is at most one term in P“^ (<pj for which 
we do not have access to pointwise estimates, or we have |p| +1 > AT-(9/2+1), 
in wich case fc<AT-(9/2 + 2) and we can bound pointwise P^ (<p) as 

If’dpfi/OI 5 (l + log(l+ t)) d . 

Pdri^\dx‘Z(<p)\zP (/)| can then be estimated as before, using Holder in¬ 
equality, pointwise estimates on p(|Z^(/)|) and the estimate 04) . 

2. Estimates on I 2 . Using 130) again, we have 

i^ZZ Z \p p d ‘m\d x >z r m\YHf)\\Y a m 

d=0 i=l |yl<|/5|, 

\V\S\P\ 

where P^(<p) are multilinear forms of degree d and signature k < |/3| - 1 sat¬ 
isfying k + |y| + |p| < \p\ + 1. Since \/3\ < 9/2 + 3, we have \f5\ - 1 < 9/2 + 2 < 
N - (9/2 + 2). Thus, using the bootstrap assumption [42] all the terms in the 
decomposition of each of the PfZ ; (<!>) can be bounded pointwise. Thus, we 
have 


u 2 i<(i+iog(i+t)) l/il+ 1 i; £ \d X iZ^m\Y\f)\\Y a m. 

i=i w<w\, 

|J?I<I/3|, 

lrl+l»?lsl/3|+i 
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Now since |y[ < \/3\ < 9/2 + 3 in the above sum, we have [y| < N - (9/2 + 1) 
since N > 14 and thus, we can bound [d_ t , Z r (0) | pointwise using the boostrap 
assumption fTTTl . We have thus obtained 

e 1/2 

ll^2(flllll(R'j)xL2( R «) < — 1 + 5 , Y ||L < d</>)(f)L^(/)(f)|| I l (K;(XR[)) , 

U + rJ \r]\<\p\ 

for some S' > 0 . 

3. Conclusions of the proof of the lemma: 

Let 

fm = Y Y \\Y a w(t)Y\f)m LH u M 

\r]\<\p\\a\<N-l 

= £ Y Y \\YiY a (cp)(t)Y\f)(m LHUM 

\r]\<\p\ i=l |a|=JV-l 

FSit) = Y Y Y \\d x iY a {cp)(t)YHf)m LHKM 

|J)|<|^| i=l |a|=JV-l 

F — F\ + f *2 + F%. 

Combining all the ingredients above, We have obtained that, there exists some 
S' > 0 such that 

r-t £. 1/2 

F{t) < / -— F(s)ds + C a e 3l2 f. 

Jo (l + s) 1+d 

Applying Gronwall inequality and using the smallness of the initial data then 
finishes the proof. 


□ 

Using (8), we have similarly 

Proposition 6.3. Letcr > 0. Then, there exists a C a > 0 such that for any multi-indices 
a,p,\a\< TV-1, |/3| < 9/2+3 and any cp £ Jl, we have for all t £ [0, T], andalll < i < 3, 

11(1 + v 2 )Wm Y a lcp)lt)Y l) {f)m L M {RicXRv) < Cad + tfe, (46) 

11 (! + v 2 ) ftf Y t Y a {(p)(t) Y? [f) It) | | £ i + « (RxX r ii) < C a (1 + tfe, (47) 

lld + V 2 ) 2 ( 1 +S) Q xi Y a lcp){t)YP{f)m L i + s {RxxRv] < C a {l + tfe. (48) 

Proof. The proof is almost identical to the proof of the previous proposition and 
therefore left to the reader. □ 

Finally, we can get rid of the small growth provided with look at a product of the 
form Y a d x {<p)YP{f). 

Proposition 6.4. For any multi-indices a, ft, |a| < AT-1, |/1| < 9/2 + 3 and any cp £ Jl, 
we have for all t £ [0, T] and all 1 < j < 3, 
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\\Y a {d x jcpmY^f){t)\\ LHRxXRv) 


< e , 


as well as 


||(1 + v 2 )^Y\d xi y){t)YP{f){t)\\ L , +s ^ xmv) <e. 

Proof. First note that the previous arguments used in the proof of Proposition [672] 
still apply and that we may only focus on the terms leading to the t a growth in the 
proof of Prooosition l6.2l which were contained in the error term 2 - More precisely, 
the term leading to the t a growth is the first one on the right-hand side of (45), i.e. 
the term tZ a d x iZ((p). In our case, this term should be replaced by 


tZ a d x id xj Z((p), 

with | a | < N- 1. Commuting the d x and Z“, we have 

\tz a d x id xj zm<t £ \d x Z v (p\. 

|j;|<iV 

We can then repeat the previous arguments (i.e. use Holder inequality and the 
bootstrap assumption (39)), except that we have gained a 1 It factor. This gain now 
means that the resulting error will decay like II t 2 ~ a , which is integrable in t and 
therefore does not lead to any growth. The L x + „ weighted estimates can be treated 
similarly. □ 

6.6 Improving the bootstrap assumptions 

We are now in a position to improve each of the boostrap assumptions. 

6.6.1 Improving the estimates on Z a {(p ) 

Similarly to Lemma l5Zl we can improve assumption (39) to 

Lemma 6.7. For all 0 < S' < 8, there exists aCgi > 0 such that for all multi-index a 
with | a | < N, 

\\S! 2 Z a <pm L i + s’ <C 5 -e. (49) 

Applying the Gagliardo-Nirenberg inequality, we deduce that for all t e [0, T\, 

||VZ“ 0 (r)|| L < 7 ( R «) <e, 

for all 3/2 < q < 1 ■ 

Proof Using the commuted equation (32) for Z“</>, we have 

11AZ“(</>)II il+ j; <J\+Jl 


where 


\a\ |a|+l 1 , , 

aT I E Hip i e > 1 

j=l d=1 |/i|<|a| 


k 1+s '’ 
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where the Up) are multilinear forms of degree d and signature less than k satis¬ 
fying 

k < |a|, k + \p\< |a|, 

and, in view of Lemma l6~4l such that when \a\ - N the only top order terms in 
Pjp (<p) are of the form Y t Y^ Up) or d x , Y^Up) with /i = N— 1, so that we can ap¬ 
ply Prooositions l6.2l and l6.3l 

and where 

\a\ , , 

/2=E E \\p[QdB^m p (g)}\\ L i^, 

d=0\p\<\a\ 

where the Q^^(d x (p) are multilinear forms of degree d of the form l33t and signature 
less than k' satisfying 


fc'<[a:|-l, k' + \p\ < \a\. 

Following the strategy of the proof of Lemma [5721 we see that it is sufficient to prove 
L p bounds on p(/i) and pUf)- With this in mind, recall that, for i = 1,2, 


llp(/()ll L P = 



< 



|| xM Vq Ji 


L 


P , 

X, V 


for any weight function Choosing x(f), p and q as in the proof of !5.2l we only 
need to prove the v-weighted L x v bounds for J\ and fa. 


1. Estimates on J\. Since k < \a\, there can be at most one term in the decom¬ 
position of each of the for which we do not have access to pointwise esti¬ 
mates. Thus, we have 

ii 5 ^(<p)^(g)i<a + a°g( i + f)) d “ 1 E E 

\ri\+\p\<\a\<P € Tm 


Now in the above sum, either \q\ > N- (9/2 + 2), in which case /j < 9/2 + 3 
and we apply the product estimates of Proposition [673] or \q\ < N - (9/2 + 2), 
in which case we can still estimate | Y n Up) \ pointwise. In this case, we use the 
^-weighted bounds on YP (/') contain in the norm En, see [20]. In both case, 
we can absorb the f-growth thanks to the t weights in the definition of J\. 

2. Estimates on Jz- These are obtained similarly, using Prooosition l6.4l instead of 
16.31 since there is no t weight in J 2 to absorb any t growth. 

□ 


The preceding lemma improves the boostrap assumption (39]. 

Similarly to Lemma l5T3l and Corollarv l5.2l the pointwise estimates on p (| (/) |) 

can be transformed into pointwise estimates for VZ“0. 

Lemma 6.8. For all multi-index a such that\a\ < N -3 

\\yZ a <P\\L 2 (R n ) ^ 2 > 


and for any multi-index\a\ <iV-(9/2+l) and Z a eT^, we have for all t e [0, T], 

e 


|VZ“0| < 


(1 + |x| + t) 312 t 112 
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Proof. We use again the commutation formula (32l . Apart from the top order terms, 
we can estimate all the quantities on the right-hand side using the pointwise esti¬ 
mates J42) on Y p {y>) for |p| < N— (9/2 + 2), (43} on Y p {dq >) for |p| < N- (9/2 + 3) 
and Proposition l 6 .ll Once we have access to pointwise estimates, we can just follow 
the strategy of the proof of !5.3l Thus, the only difficult term are thos containing top 
order terms in Y p (cp). Those coming from the P multilinear forms are of the type 

~ t P{Y p «p)f) 

with |p| the largest integer such that |p| < N - (9/2 + 1), i.e. |p| = N - 6 . Now, since 
|p| < IV- 6 , we can consider Y^{Y p {(p)f) for \/3\ <3. Thus, we may apply Proposition 
I6.1l directlv to the product Y p (<p)f and we find 


7 P[Y p (cp)f) 


< 


e 

t l ~ a (1 + f + |x |) 3 ’ 


where we have used Proposition 16.21 to bound the norms appearing on the right- 
hand side after applications of Proposition [64] Choosing er < 1, these terms there¬ 
fore decay better than what is needed for the statement of the Lemma. 

The top order terms coming from the Q forms are of type 


p[Y p {d x i<p)f), 


with |p| < N- (9/2 + 2), i.e. |p| < N - 5. We can then proceed similarly. Since we 
have no extra t decay in front of the Q forms, it is important not to lose any t decay 
here and thus we use the improvements of Proposition l6~4l The rest of the proof is 
identical to that of Lemma [5731 and therefore omitted. □ 


The preceeding lemma improves (40) . To improve fTTI is then suffices to use 
Lemma [OH as well as the pointwise bounds on Y a Up) LI 21 . Thus, it remains only to 
improve (38), (42) and (43). 


6.6.2 Improving the global bounds 

We have 

Lemma 6.9. For all t e [0,T], 


E N ,g[f{t)] < 3/2e, 

which improves (38) . 

Proof. Using the commutation formula (30) . we have 

[^,F“](/) = '“£ 1 f ^ P“’jj<p)d x iZV{<f))Y p (fY 

d=0 /=ll7l<|a| 

\P\<\a\ 

where the (</)) are multilinear forms of degree d and signature less than k such 
that fc<|a|-l<iV-l and k+\y\ + \f \ < |cr[ + 1 <AT+1. When k< N - (9/2 + 2), we 
can bound P“ p g Up) by a polynomial power of 1 + log(l +1) and we can treat the other 
terms as in Section [534] Otherwise 1c > N - (9/2 + 2) and thus \f>\ <9/2 + 3 and we 
can use the estimates on products of Proposition l6.2l As in Section l53l it follows that 
all terms are integrable in L 1 (R, ; /J (!” x R'J)) and the lemma follows from Lemma 
[ 22)1 □ 
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6.6.3 Improving the pointwise bounds on Y a (ip) and Y a (d x , <pj 

Finally, we conclude the proof of the 3d case by improving l4~2l and 1431 . 

Lemma 6.10. ForallipeJi andforallte [0,7’], 

1. for all multi-index a with\a\ <N-{ 9/2 + 2), 

|F“((p)|<c(l+log(l + t)), (50) 


2. for all multi-index a with |a| < N— (9/2 + 3) and all 1 <i<n, 

| Y a {d x up) \ <e. (51) 


Proof. Using the commutation formula l30l . we have 

[7>,F“](<p) = l “£ 1 f ^ 

d =0 i= 1 |y|<|a|, 
l/3|<|a| 


where the P^p (<P) are multilinear forms of degree d and signature less than k such 
that k < |a| —1 < AT-1 and k+|y| + |/l| < |a| + l < Af+1. Given the range of the indices, 
we can estimates all terms on the right-hand side pointwise and find that 


|[r 0 ,F“]((p)|< 


(l + log(l + t))'“ l+2 c 
t 2 


which is integrable in t. On the other hand, we have Y a T f i, (<p) - Y a fV Z[<p). Using 
Lemma fbBl it follows that 


|a| 

Y a T^cp) = tZ a VZ(l)+ £ PUcp)zNzf, 
d =1 

where the second term is integrable in t, while the first term satisfied only the weak 
bound, for all 

\tZ a VZd>\< 

r 1 + t 

Since any solution to ff g) - h with 0 initial data satisfies, for all t > 0, and for all 
x,ve IR” x 05”, 

| g\(t,x,v)< / sup \h{s,x,v)\ds, 

Jo (jc'.y'JeRjxRg 

l50l follows. FTil can be obtained simarly, using that 

\tZ a SIZd x np\ < 

since d x ,<p= ~ t [td x i) {([>) with td x i ef. □ 
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